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Analytical Results in Magnetic Saturation 
Amalityczne rezultaty dotyczące wysycenia magnetycznego 


AHaJlkTUTECKME peLlIeHMA NPOÓJEMEI MATHMTHOTO HACHIUJCHMA 


1. INTRODUCTION 

Saturation of magnetization is of interest in many applicattons and in 
particular in application to magnetic recording heads. Current work on 
calculation of this nonlinear effect uses the direct numerical integration 
of the nonlinear equations [1] The purpose of this paper is to explore 
analytical methods in the hope that they will be useful in uderstanding 
this and other similar problems. 

We analyze the solution, in two dimensions, to the equation V-B=0, 
where B = +1V9, and u is a function of H = vel. It is therefore also 
applicable to the case of dielectric saturation, where B is replaced by 
E and u by €, as well as steady two dimensional flows of an Ideal 
compressible fluid where u is replaced by p, the fluid density. We shall 
however focus on the case of magnetic saturation effects with the 
eventual possible application to recording heads. 
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Since the Linear problem is solved by the real or imaginary part of 
any analytic function, it is naturel to ask if any such harmonic solutions 
are possible for the nonlinear case. The answer to that is by no means 
obvious. We show however that such nontrivial harmonic solutions do 
not exist. Although there is no reason to expect such solutions to be 
possible, it is perhaps surprising that even the conjugate potential 
cannot be harmonic for any but the trivial case. 

We explore the possibility of solutions which are parallel to 
harmonic functions. Since even in the nonlinear problem the potential 
and field Lines (or conjugate potential) are orthogonal, it might be 
possible to have nonharmonic solutions which are parallel to harmonic 
functions. We find that this is possible but for only two special 
solutions. Thus we find two exact nonlinear solutions and suggest an 
approximate method based on these. 

Exploring further a complex function formulation, we use the 
complex potential of the Linear problem to conformally map the x-y 
plane to the a-B plane. The equipotentials of the Linear problem are a 
= constant, and the field Lines are g = constant and are thus coordinate 
Lines. Saturation effects will distort these straight Lines. If we assume 
that the distortions are such that the nonlinear equipotentials in the 
a-B plane have small slopes, we can find a Linear equation for a 
certain “v-potential". Once this equation is solved for v, the first 
denvatives of the potential $ are related to v through an algebraic 
nonlinear equation which can easily be solved. 

In section 2 we set up the equations for the model, and in section 
3 introduce the conjugate potential. Section 4 interprets the meaning 
of possible harmonic solutions. Formulation in terms of the complex 
variable is developed in section 5 where the Linear and nonlinear cases 
are contrasted. In section 6 we show that the conjugate potential 
cannot be harmonic. Section 7 develops the exact solutions which are 
paratlel to harmonic functions and concludes with the deduction that 
the potential cannot be harmonic. Section 8 suggests an approximation 
method using weak dependencies. Another approximation method is 
introduced in section 9 where Liner approximate equations for the 
v-potential are derived. 
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2. THE MODEL 

The magnetic permeability H, of ferromagnetic matenals is 
generally not a constant but depends on the field strength H. This 
dependence is such that the nonlinear effect is called “saturation”. 
Thus as the field strength increases, u decreases. A simple, commonly 
used model for this dependence is Monson's model [2] which is given 
by: H = Holl + x/(1 +H/H,}J, where po, x, and H, are constants. 

Typically for magnetic cores of recording heads x ~ 10°, and Hel 
Oerstead. Saturation is strongest near the comers of the recording 
heads where H is Largest and according to numencal calculations |1| 
reaches orders of 500H.. Therefore even at and above saturation 
fields the x term in is 10 times es big as the 1. Thus neglecting the 1 
term and using H, as the unit of the field we can wnte approximately: H 
= Hox/(1 + H). These approximations are not necessary, they are made 
for simplicity. The methods developed here will also apply to other 
models of ju(H]. 

The equations to solve inside the magnetic medium are therefore: 


H=-Vó and V-[UV$]=0or V-fV$/(1+ HY -0 (1 
which is the same as V-[1V$]= 0 with u replaced by V(1 + H). 


3. REAL FUNCTION FORMULATION 
We confine ourselves to two dimensions. Define V* = (-8/2y, 8/8x]. 
This operator is “orthogonal” to the normal "unstarred" gradient. That 
is VV” = V*V = 0, end is equivalent to the operator -kxV, where k is 
the unit vector perpendicular to the x-y plane. 
The following properties of this operator are evident: 


V**=-U ; (U*y]Vy =-(Uy)-V*v=0; V*-U*=VŻ; |U*y|=|Uv|. 
Therefore it follows that when V-B = 0, we can find a scaler $ 
such that B=-V"y. (The negative sign is not significant, but is picked 


for Later convenience.} Using this we get from (1): 


HV =-V*y (2) 
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It follows from {2} that the curves d = constant and y = constant 
are everywhere orthogonal. This in tum means that y is constant 
along field Lines, which are orthogonal to the equipotentials. In 
analogy with the Linear problem we call 4 the "Conjugate Potential”. 


4. PHYSICAL REQUIREMENTS FOR HARMONIC SOLUTIONS 
Consider the general case where y is not harmonic. We obtain a 
physical condition for a harmonic potentiai by evaluating V7y. Operate 
with V*- on both sides of equation (2] to get: 


V*[HV$) VY =-(V*H)-Vó LVV p/p (VYE 


where we have used [2} to replace Vo with - V*y/u. We have 
also substituted © =tnlu]. Therefore 


UŻy =-(V£)-V4 (3) 


Similarly using the starred version of equation [2] and operating 
with V- on both sides, then making similar replacements we get 


726 = (VL Vb (4) 


Consider the condition necessary for harmonic solutions. The RHS 
of equation (3) will be zero only if € = constant fie. i = constant) or VE 
and Vy are orthogonal. Since we presume that p= constant, the only 
way to have a harmonic solution for w is to have £ = constant and y = 
constant be orthogenal. Since ¢ = &(H), it follows that this condition 
requires that the magnitude of the H field change onty along the field 


Lines and not perpendicularly. 
SimiLarLy from equation (4) we see that if $ is to be harmonic, then 


the H field magnitude should only change in the perpendicular direction 
to the field Lines and not parallel to them. 


It will be shown that neither of these cases is possible except for 
the trivial case of a uniform field. 
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5. FORMULATION USING FUNCTIONS OF A COMPLEX VARIABLE. 

Using Z = x +Iy, with z* = x - iy, and treating z and 2” as formally 
independent variables, we can deduce that 3/82" = 3{8/x + i8/8y]. 
Therefore since u,$ and # are all real functions, the two components 
of equation (2) can be written as the real and imaginary parts of 


186/82" + i 3y/8z*=0 (5) 


If u = constant, then equation (5) can be wntten more simply as 
alué +iy)/Az* = 0, which implies that 


Ho + ip = fiz) (6) 
where f(z} is an analytic function. 


This is the Linear case in which any anatytic function of z is a 
solution of the equations, where the real part is id and the imaginary 
part is +, the conjugate potential. 

When u is not constant, however, we can always find a function n 
such that 


„34/0z2* = On/6z* =- 104/02". (7) 
From which it follows that 
Aln + i¥)/8z* = 0 and thus n + ty = F(z) (8) 


Note however that whereas if u is constant, q = Hb, and hence is 
real, when p is not constant, n is not reat. In fact all we can say at 
this point is that it need not be real. The fact that it cannot be real 
except for the trivial special case of a uniform field, follows from the 
following. 

Suppose 1 is real then it is harmonic because by (8) it is the real 
part of F(z] and 4 is its imaginary part. Transform (7) from the pair z,z* 
to the pair F, F* as independent vanables to get 


La$/2F* = an/2F* 
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or in reat terms 
Hêé/an-1, and 18%/0y = 0; (9) 


Hence $ = ln) and from the first of equations (9) we conclude that 
p= pln). But p= (1+ |F'|o'r’, where pnme denotes differentiation. 
Therefore [F'| = fn]. Taking this to its Logical conclusion we deduce 
that since inff) is harmonic 8*lnff]/an? = 0, and so f = |F'| = fe. This 
Leads to the result that F' = Ae, and so F' must be a constant, hence 
the trivial solution. This proves that n cannot be real for nontrivial 
solutions. 

Proceeding with the general formulation we note that (8) relates 
the nonlinear problem to an analytic function F(z). However only the 
real and imaginary parts of F(z) are harmonic. Since n and ẹ are not 
the real and imaginary parts of F, they are, in general, neither 
orthogonal nor harmonic. 

Companng the linear and nonlinear cases we note that in the 
linear case given an analytic function we immediately have the 
potentials. In the nonlinear case on the other hand, given the analytic 
function F{z) we still do not know the potentials. 


6. w NOT HARMONIC. 

Consider the possibility of the special case Vy = 0, ie. y is 
harmonic. We show that n can be made reat. 

Since by (8) n= F(z) - ty It is also harmonic. We can always add an 
arbitrary analytic function of z to n and not change the determining 
equation [7]. Therefore without Loss of generality, by adding anatytic 
functions to n we can arrange to make it real at which point n = RelF[z]] 
and y =Im{F{z}} Since we have shown however, that n cannot be real 
for nontnvial solutions, it follows that y cannot be harmonic. 


7. EXACT TIONS 
Equations 3 and 4 indicate a great deal of symmetry between the 
‘potential $, and the conjugate potential w. It is not surprising 
therefore that $ cannot be harmonic either. Rather than prove this at 
this time we shall derive two exact solutions, which would include 
harmonic ones if they existed. 
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Since we have ruled out the possibility of harmonic solutions for y 
of the nonlinear problem, it is natural to ask if there are any solutions 
which are functions of the harmonic Linear ones. 

Consider the fact that if a,B are the real and imaginary parts of 
f(z] and so harmonic, nonlinear functions of a and p (say fila) and 
f‚(g)l, are not harmonic. We Look for exact solutions to the nonlinear 
potential $ in the form of a function of the real part of f(z). 

These exact solutions may be useful as starting points for more 
general solutions. In particular we suggest how the use of “slowly 
varying coefficients” can generalize these exact solutions to 
approximate, but more general, solutions where the equipotentials a = 
constant Locally approximate the actual equipotentials. 

Let a and B be two harmonic functions such that F(z] = a + IB. We 
pick œ and 8 such that the curves a = constant and $ = constant can be 
made paratlel, and therefore 6 = constant and w = constant are also 
parallel. It follows that $ = bla} and y = w(gl. Substitute this into 
equation (2) and obtain: 


Hó'Va =- w'V*B, where $' = dó/da etc. 
Writing this in complex notation we have 

Hip'/y'ldax/8z* + i8g/8z* =O (10] 
But since ala + ip)/dz" =O, it follows that 

plo'/y')=1 (11) 
Substituting p = {1 + HY! we get: 

H/o' = Vy! - 7/9" (12) 


Substituting for H= |V6| = |$'Va|= |6'F'(z]| and Letting € = sign(4') ( 
= | for 60 and - for $'<0), we obtain 


|Vale=Vy'- 1/6" = |F(z)|e (13) 
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Since € is either +1 or -1, if we solve equation 13 with e= 1, then €= -1 
case is simply y and $ interchanged. This presumes that € does not 
change sign inside our region. Clearly € can only change sign when d' 
vanishes. But $' can only vanish at a singularity of F. We assume that 
this can only happen at isolated points in our region, and therefore no 
transition of signs is possible. Consider therefore €= 1. 

We examine the conditions under which this equation has sotutions. 
Take the Logarithms of both sides of equation 13: 


Ln|F'(z)| = LnlA(a] + Bie) (14) 


where Ala) = -Vé' and Bip) = 1/4'. Since F'{z) is analytic, Ln|F'| is 
harmonic. Therefore V? Ln|F'| = 0= V@Ln(A + B). Therefore 


(82/8a7 + 867/882] Ln(A + B}=0 (15) 


Equation (15) is clearly both a necessary and sufficient condition 
for (13) to be satisfied. Evaluating the derivatives we have: 


A'A + B'B + AB" + BA" =(A'F + (BE (16) 


If (16) is to be true then 82/8aa8 of the equation is also true. 
Therefore a necessary [though not sufficient) condition for (13) is: 


A'B" + BA"=0 (17) 
Separating variables we get: A"/A' =- B"/B' =A = constant. 
We have 3 possible cases: 
A=0 
A> 0, typify with A =1 
A <0, typify with A =-1 This case however is simply a<->8 
exchanged and so needs no separate consideration. 


(Different values of A amount to different scalings of a and 8] 


Consider A = 0, then A" = B™=0. 
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Rap + aa + 8207. 
B= Do * DB + Dap”. (18) 
This solution has to be substituted back into equation (16) in order 


to get the necessary and sufficient condition for a solution. When this 
is done we obtain 


82=bz; A; + b? = 4[ag * bolaz. (19) 
Considenng the case A = 1 we get B” =-8' and A" =A’, therefore 
A= + 6e” + 28%; B= bo + bicosB + Dosing. (20) 


Substitution of these back into equation 16 gives the necessary 
and sufficient condition: 


Do =-8o; Di? + bo? = 46,82. (21) 
We now solve for $ and y, and then for Ffz]. 
Consider first the A = D case. We can simplify the algebra 


considerably by noting that if we shift the ongin of [a,B] we can 
arrange to have the quadratics in (18) without the Linear terms: 


A= Gg + aa. Where a, = a + ay/[262]; ag’ = (48962 - a;7]/1462] 
similarty 
B =bp' + bzpy?. Where Bj =B + by/[2b2); Do’ = (4bobz - b;2)/14b2] 
Therefore define F(z) = a, + if, = F(z} + c, where c = - (a/[2a2]* 
b,/(2bz]) Then |F'y| = elan? + py?) = 87 | Fi |? or Fy’ = agf}, where y is 


a real constant. Integrating, this gives 


F=- SYBo/(Z - Zp). (22) 
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Without Loss of generality, picking the ongin at zo and choosing our 
axes appropriately (7 = 0) we get F,=-82/z. where a> is real. Now a, = 
-87Xx/(X? + y?), and p; = -82y /L< + y*), and so the equipotentials ay = 
constant and the field Lines By = constant are circles. The field is given 


by 

H=V$=(Vaq]$'=- Va,/A= -Voa/lag' é 8207) (23) 

and the $ and y potentials by 

$ = 9 - foan Bo * Bia + 820?) = po * Ln - a,]/ (a - ap 

= bo - 2k tan (k/(2eza + ay) 

Y= bo * dB/l Do + Di + D2B7]= Yo + k"Lni(B - B, (B- BJ) 

where 2a, = -84/67 + vk/82, 2B, = Dy/87 + k/a3, 

and k = /l48982 - a;7). 

This solution requires & = 0 and represents field Lines with the 
same shape as those of a Line dipole in the Linear case. 

For the special case 4 = 0, equation 19 requires a; = 0 = by also and 


the result is = 9 + œ/8 and y= bo + B/Do, Le. the harmonic solution. 
Substituting this into the equation for F(z) we find in that case 


F(z) = [ag + bolz + C. 
or the tnvial case of a constant field. 
The case A = 1 is similarly simple. Rewrite A and B from equations 
(20) as: © À 
A= ag + 2/lejgzlcoshla + 7) 


and B =- 69 + 2yla;G>lcoslB + 74), (24) 


where tanhy = (a, - 6>]/ (a, + as] and tany; =-Dp/D). 
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This gives the result: 

[F"|= 2 laraz]l cosh(a,) + coslpy) } (25) 
where we have substituted 

Q,=0 + y end B =B + 7, ANdF,=F +(7 + 17) 
Equation 24 may also be wntten as 

JF, |= 2V(a;82] | cosh(3F,) |? (26) 


This is solved by: F', = 2ylajazle"cosh"[1F,], where © is a real 
constant. The solution to this in tum gives: 


Fi“ Ind + zR]- Ln(l- zR) (27] 


where we have replaced /lajazle® by R, allowing R to be complex. 
z may be replaced by Z - zo, where Zp is an arbitrary complex 
constant. 

The potentials are given bu: 


$ = ġo - ky 'Lnffe**7 - a,]/[e**7 - a_}} 
where a, *-ag/r + ky/r; r=2yla;82); k,=Vla3ż-r7] 


The equipotentials for this solution are also circles. For the 
special case of 87 = 0, equation 19 again requires 4, = O = b, which gives 
the trivial result of F = constant, and thus zero fietd. 


We have deduced all the solutions to the nonlinear problem 
whose potential can be written as a function of the real part of an 
analytic function F. This includes the possibility of harmonic solutions 
for the potential $. 


Consider now the possibility of a harmonic solution in the Light of 
our results. Since = ola) then V76 = 6'V2a + [VaF$" = (Valo" = 0 
for a harmonic solution. Therefore " = O is the only possibility. The 
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only such solutions we have found (both A = 0 and A = 1) are the trivial 
constant field solutions. 

In conclusion therefore we have shown that the only harmonic 
solution to the nontinear problem is the trivial one of a constant field. 

8. THE ADIABATIC APPROXIMATION. 

Suppose we now relax the condition that $ = constant be exactly 
parallel to a = constant everywhere, and replace it with the condition 
that it be approximately parallel everywhere.. This implies that $ = 
o{a,B) where the g dependence is weak in some sense. This implies 
that the denvatives (using A$/3a = 4, etc) |o.]> [dp]. 

For a graphical picture of this condition imagine the curves $ = 
constant drawn in the a, B plane. The condition then requires the 
curves 4 = constant are to have small stopes, l.e. to be approximately 
parallel to the g axis. Similarly we can require y = v[B,ax] where the a 
dependence is weak. Since we must still require orthogonality between 
the $ = constant and y = constant curves, it can easily be shown that 
our derivation from equation 10 through 27 is unchanged, provided we 
drop all second order derivatives with respect to the weak 
dependencies. This allows us to proceed with more general problems. 


2. DERIVATION OF APPROXIMATE EQUATIONS. 

Suppose ‘ve use the Linear solution as the starting point for the 
nonlinear one. We can then use the potentials [both the real potential 
and the conjugate potential) as our "grid" or curvilinear set of 
coordinates for the problem. This will allow us to assume that the 
actual potential curves have small sLopes in terms of our coordinates. 
This in turn, as we shall see below, allows us to derive Linear 
equations for a “v-potential”, even though the equivalent equation for 
the potential $ is nonlinear. 3 

It is important to realize that a curve with small slopes can 
deviate considerably from a straight Line. Thus in this approximation 
the nontineer potential can differ considerably from the Linear one. 

We start with equation (5). Using w = F(z) = œ + ip [and with an 
obvious notation, the implied relation w* = F*(z) = F(z*)} we transform 
equation (5) from the independent variable pair z,z* to the pair w,w*. 


}0o/dw* =10y/0W* = Ony/dW*. (28) 
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In (28) we substitute n= a + iv, lor the equivalent 4 = B - vl and 
equate the real and imaginary parts separately to obtain 


Hobe =1- Vp, Hóg = Vee; (29) 


where the subscripts indicate partial denvatives with respect to 
the independent variables a and B. To comptete the set we add to (29) 
the expression for H: 


use |F'|(62 + 63211, where F' = dF/dz. 


The following procedure will yield the exact equation for v: The 
pair of equations (29) is solved for.$ and ds in terms of vy and vy. 
Then the "Consistency Condition" x, = ds. IS applied giving the 
equation for v. Of course the alternative procedure of applying the 
Consistency Condition Veeg = Vp. to equations (29) will yield the original 
equation for b with a,8 as independent variables. At this point there is 
nothing to choose between these two approaches, and thus no 
simplification apparent. The following discussion will make clear how a 
judicious choice of F(z) will allow us to make reasonable 
approximations and dictate the use of the equation for v rather than 9. 

Let us choose F(z) to be the complex potential of the Linear 
problem. Thus the Linear solution, indicated by the suffix "L", is given 
by é = @ = mM and h = B, with n =0. In the ap plane the 
equipotentials are parallel to the B axis, and the field Lines are 
parallel to the æ axis. Since the nonlinear solution will in most of the z 
plene follow closely the Linear one, except in regions Like comers, 
where saturation is important, it is reasonable to suppose that in the 
aß plane the approximate nonlinear solution has a small v. From 
equations (29), and the fact that dj, = I and ds = O, it follows that 
|| || is a reasonable approximation. Therefore as a first 
approximation we can deduce a Lineanzed equation for v, or an 
equivalent one for $. However that the equivalent approximate 
equation for dis ncnlinear and rather involved, whereas the equation 
for v is Linear and considerably simpler. 
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Proceeding as outlined above, and neglecting nonlinear terms in v, 
we obtain the following equation: 


Veel!” f) + Veel” 2f] + Vafa + 2 Vels” fa. (30) 


where f = |F'|. If we further assume that f is smalt (which Limits 
the approximation to weak nontineanties) we get the much simpler 
equation: 


Vera + Vag sie fp. (31). 


Equation {31} is particularly easy to solve, since the solution can 
be written in terms of the Green's function Glx‚y.x'‚y') for the Linear 
problem. To exhibit this we re-write equation [31] in the z-plane as 


UŻy = - fg/f? = 8['1)/8B 
which has the formal solution 
vix,y} = [ [G6lxy,x'„yJ8fF10e„yJ]/Opdx'dy' 


The Linear solution is assumed to satisfy the correct boundary 
conditions and so we have used the fact that v = o on the boundary. 

Since the Green's function G is easily determined from analytic 
function theory (see for example Hildebrand 1976), this part of the 
problem is standard. 

Once we have the solution for v, the equations relating $,, and be 
to Vv, and vs [equations 29] are algebraic, though nonlinear, and can be 
easily solved. a 


10. CONCLUSION. 
We have shown that there are no nontrivial harmonic solutions to 
the nonlinear problem. Exact solutions heve been found and 
suggested as the starting point for a general approximation method. 


A new potential function [v} is defined which approximately 
satisfies a Linear equation (equation 30]. 
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In a future publication these equations will be solved for some 
realistic geometnes. 

The methods developed in this paper should also be applicable to 
other physical nonlinear problems. Electrostatic saturation is an 
obvious case. The two dimensional ideal compressible fuid flow 
equations though somewhat different in detail, should be amenable to 
a simmilar treatement. 


REFERENCES 


1 Bertram Niel H. and Steele Charles W. IEEE Trans. Mag. 1976, 
MAG-12, 702-206 
2. Monson J. E. IEEE Trans. Mag. MAG-8, 536-538. 


| 
STRESZCZENIE 


Metodami analitycznymi badano problem nasycenia magnetyzacji 
w dwu wymiarach używając modelu Morisona dla p Zdefiniowano po- 
tencjał magnetyczny i potencjał sprzężony i pokazano, że nie ist- 
nieją dla nich nietrywialne rozwiązania harmoniczne. Yykazano, że 
istnieją tylko dwa dokładne rozwiązania o liniach ekwipotencjal- 
nych równoległych do potencjałów harmonicznych. Rozriazania te 
mogą stanowić punkt wyjścia dla przybliżonych ogólniejszych roz- 
wiązań, 

Zaproponowano także metodę przybliżeń wykorzystującą funk- 
cje zespolone i prowadzącą do liniowego równania dla "potencja- 
łów VY nawet rtedy gdy odpowiednie równania dla rzeczywistego 
potencjału są nieliniowe. 


PESOUE 


C MOMOMBN AHANMTMYeECKOrO METOĄA, NPA UCHOLb3OBAHAM MOJENN 
KopucoHa ANA m, DPaCCMaTpPMBaNach MpoOleMa HACHIleHMA MATHETU- 
3aNKU B ZBYX pasMepax. OnpSedeleHH MAFHUTHHi M CONPAREHHHA NO- 
TEHUHAJIH M NOK43AH0, UTO ALA HEX HE CYWECT2y0T HETPUBKAJHBHHE 
TapMOHUUECKHE peueHuR. JokastBaerca, YTO CYHECTBJKT BCETO ZBa 
TOVRNE PeweHMA C ZKBENOTEHLAARBHHMA AMHAAMK NAPANIENBHAMA rap- 
MOHMYCCKAM NOTEHHMA1AM. GTM PELEHUA MOTYT ÓNTB HCXORHHMH ANA 
NpPMÓJKEEHKNK, COTES OÓiAX pemeHXń, 

Ipezroxken TORE METOĄ HPAÓJAXSHKH UCHONBZYDNKK KOMNIIEKCHNIE 
Qyskunu, Bezjuuit K NMHEHHOMY YpaBAEHUK ana "V -NOTEHNMANOB" 
make TOTZA, KOTZE COOTBETCTBYDNAE YDABHEHKA ANA AZeACTBUTeNBHOTO 
norengana - H6NAHGHHNE, 
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